Accepted in ApJ 

Preprint typeset using I^T^jX style cmulateapj v. 10/09/06 



THE "TRUE" COLUMN DENSITY DISTRIBUTION IN STAR-FORMING MOLECULAR CLOUDS 

Alyssa A. Goodman 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, MS-42, Cambridge, MA 02138, USA 

Jaime E. Pineda 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, MS-10, Cambridge, MA 02138, USA 



OO 
O 

o 

(N 
•*-> 
o 

O 



43 

6 



> 

m 

\6 
o 
oo 
o 



X 



Scott L. Schnee 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, MS-10, Cambridge, MA 02138 and 
Division of Physics, Mathematics and Astronomy, California Institute of Technology, 770 South Wilson Avenue, Pasadena, CA 91125, 

USA 
Accepted in ApJ 

ABSTRACT 

We use the COMPLETE Survey's observations of the Perseus star-forming region to assess and in- 
tercompare three methods for measuring column density in molecular clouds: near-infrared extinction 
mapping; thermal emission mapping in the far-IR; and mapping the intensity of CO isotopologues. 
Overall, the structures shown by all three tracers are morphologically similar, but important dif- 
ferences exist amongst the tracers. We find that the dust-based measures (near-IR extinction and 
thermal emission) give similar, log-normal, distributions for the full (~ 20 pc-scale) Perseus region, 
once careful calibration corrections are made. We also compare dust- and gas-based column density 
distributions for physically-meaningful sub-regions of Perseus, and we find significant variations in 
the distributions for those (smaller, ~few pc-scale) regions. Even though we have used 12 CO data to 
estimate excitation temperatures, and we have corrected for opacity, the 13 CO maps seem unable to 
give column distributions that consistently resemble those from dust measures. We have edited out 
the effects of the shell around the B-star HD 278942 from the column-density distribution compar- 
isons. In that shell's interior and in the parts where it overlaps the molecular cloud, there appears 
to be a dearth of 13 CO, which is likely due either to 13 CO not yet having had time to form in this 
young structure, and/or destruction of 13 CO in the molecular cloud by the HD 278942's wind and/or 
radiation. We conclude that the use of either dust or gas measures of column density without extreme 
attention to calibration (e.g. of thermal emission zero-levels) and artifacts (e.g. the shell) is more 
perilous than even experts might normally admit. And, the use of 13 CO data to trace total column 
density in detail, even after proper calibration, is unavoidably limited in utility due to threshold, 
depletion, and opacity effects. If one's main aim is to map column density (rather than temperature 
or kinematics), then dust extinction seems the best probe, up to a limiting extinction caused by a 
dearth of sufficient background sources. Linear fits amongst all three tracers' estimates of column 
density are given, allowing us to quantify the inherent uncertainties in using one tracer, in comparison 
with the others. 

Subject headings: dust, extinction — ISM: abundances — ISM:molecules — ISM:individual (Perseus 
molecular complex) 



1. INTRODUCTION 

In this paper, it is our goal to use data from the COM- 
PLETE Survey of Star-Forming Regions 1 to assess and 
intercompare three methods for measuring column den- 
sity in molecular clouds: near- infrared extinction map- 
ping; thermal emission mapping in the far-IR; and map- 
ping the intensity of CO isotopologues. We wish we 
could offer a snapshot of volume density, as numerical 
modelers can, but volume density cannot be mapped 
from our vantage point on Earth without making very 
model-dependent assumptions. We discuss the position- 
position- velocity distribution of material in other papers 
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(e.g. iRosolowskv et all 120081: iGoodman et al. I l2008h : 
here we focus solely on the column density distribution. 
Our aim is to offer insight into the biases of particular 
techniques, and to provide the best estimates to date of 
"true" column density distributions in molecular clouds 
suitable for comparison with current and future numeri- 
cal simulations. 

We focus our discussion on the Perseus molecular 
cloud complex, which covers nearly 10 square degrees 
on the sky ( ~ 200 pc 2 at 250 pc; lEnoch et all 120061 : 
iHirota et al. 2008). A companion paper to this one 
([Pineda et al.ll2008T) includes an historical perspective on 
the use of dust extinction and emission and of CO iso- 
topologue mapping to study column density in molecular 
clouds, and it offers an in-depth look at abundance vari- 
ations in the cloud. 

2. DATA: THREE MAPS OF PERSEUS 
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In the past, three principal methods have been used to 
chart column density in molecular clouds: (1) extinc- 
tion mapping, using either star coun t ing (e.g. iBarnard 




ICernicharo &: Bachillerl Il984t ICernicharo et al 
color-excess measurements (e.g. lLada et al 



or 



199J); (2) dust-emiss ion mapping, at f ar-IR through 
mm-wavelengths (e.g. lSchlegel et alll998l ); and (3) map- 
ping integrated intensity of molecular-line emis- 
sion, usually 12 CO or 1 3 CO on large (> 1 pc) scales 
(e.g. iPadoan et alj [l999h . As noted above, there is a 
long history of the applic ation of these techni ques, and 
we analyze that history in iPineda et al.l (|2008h . 

All three of these methods have now been applied, as 
part of COMPLETE, to the large (~ 10 deg 2 ) swath 
of sky in Perseus 2 which has also been surveyed by the 
Spitzer Space Tele scope under the "Cores-to-Disks" (c2d; 
lEvans et al.ll2003l) Le gacy Program 3 . The Perseus maps 
u pon which the analys is in this paper rests can be found 
in ISchnee et all (|2005l ) and lRidge et~aTl (|2006al|b1 ). 

2.1. Extinction Mapping 

The exti n ction map used here is presented in 
iRidee et all (|2006af ). To create the map, COMPLETE 
collaborators Joao Alves and Marco Lombardi, applied 
their "NICER" (Near-Infrared Color Excess method 
Revisited) method to 2MASS near-infrared maps of 
Perseus 4 . 

As described in lLombardi k. Alves! (|2001h . the NICER 
algorithm calculates extinction values based on the dif- 
ference between the observed average near-IR color 
of stars within a sampling box to the stars' intrin- 
sic average color. As the number of stars in a sam- 
pling box is reduced, the uncertainty in the extinc- 
tion calculated by NICER rises. The details of the 
tradeoffs between high-resolution/high-uncertainty and 
low-resolution/low-unce rtainty mapping with N ICER is 
discussed in detail by lLombardi et alj (|2006f ). The 
NICER/2MASS map used here is made with 5' resolu- 
tion on a 2.5' grid, a common resolution to which all data 
presented in this paper have been smoothed. As shown 
in Figure [TJi, we can measure extinctions from Ay — 
to 10 mag with typical reliability better than 0.25 mag 
(Figure [TJd) , but it remains true that the highest uncer- 
tainties (~ 0.4 mag) are in the highest extinction regions 
{Ay > 10 mag). 

Even though 5' resolution sounds coarse (especially 
compared with the 46" intinsic resolution of the COM- 
PLETE spectral- line maps) , it is important to appreciate 
that earlier, optically-based, extinction mapping meth- 
ods could never come close to this kind of resolution in 
high column-density regions. In the U-band, for exam- 
ple, even modern star-counting techniques typically fail 
(as zero stars are present in a counting box) at Ay ^ 6 
mag (e.g. lCambresyll999[ ). The material with 5^Ay^l0 
mag is exactly the material, when observed with the kind 
of few-tenths-of-a-pc resolution we have here, that ap- 
pears most actively engaged in forming stars in Perseus 

2 Ophiuchus and Serpens have also been observed similarly in 
the COMPLETE and c2d Surveys, but this short paper focuses on 
the most studied maps to-date, which are of Perseus 

li http: //peggysue . as .utexas . edu/SIRTF/ 

4 2MASS data used are from the "Two 

Micron All-Sky Survey. Final Release" (See 

http://www.ipac.caltech.edu/2mass/releases/allsky/doc/cxplsup.html ) 



(|Kirk et al.ll2OO6Hj0rgensen et al.ll2006f ). and it is nearly 
completely opaque in optical wide-field surveys. 

We choose, at this point in time, to use the 2MASS 
data for our extinction maps, rather than the Spitzer 
IRAC c2d data for two reasons. First, the 2MASS data 
cover the full region available in COMPLETE 's molecu- 
lar line maps, while the IRAC data only span "most" of 
the region. Second, and more importantly, work is still 
ongoing (Huard, private communication) on the calibra- 
tion of IRAC-based extinction maps and on the deter- 
mination an agreed-upon mid-IR extinction law. In the 
future, the best extinction maps will likely be constructed 
from photometric observations spanning the NIR and 
MIR range. Since the absorption by dust becomes pro- 
gressively less severe at longer wavelengths, the MIR data 
will be especially useful for mapping out the structure of 
very high column density gas. 

2.2. Far-Infrared Emission Mapping 

The far infrared data used here are taken from the IRIS 
Improved Reprocessing of the IRA S Survey) database 
Miville-Deschenes fc Lagachell2005t ). IRIS represents a 
substantial improvement over the earlier ISSA release of 
the IRAS survey, in that it has better zodiacal light sub- 
traction, destriping, calibration and zero level determi- 
nations. The importance of finding and using the proper 
zero level when deriving column density maps where 
flu xes are not much higher than the noise is explained 
in lArce fc Goodmanl ((2001) . For IRAS-based observa- 
tions of clouds like Perseus, the FIR signal is typically 
only much greater than the noise above Ay ~ 2 mag. 5 

Using the IRIS flux maps at 60 and 100 fim, we 
derive dust temp erature and F IR optical depth maps, 
as explained in ISchnee et alj (|2005| ). We use the 
2MASS/NICER column density map described in §2.1 
to calibrate the conversion between FIR optical depth 
and NIR extinction. Thus, we fix the overall column 
density scale with absorption measurements, and then 
adjust the dust emission properties, given derived color 
temperatures, to minimize differences beetween absorp- 
tion and emission probes of dust column. (Note that 
this procedure assures that the NIR (extinction) vs. 
FIR (emission) column desnity plots shown in Figure [3] 
have slope and intercept close to one and zero, respec- 
tively.) By comparing Figure [TJ;, which shows the re- 
calibrated, IRAS/IRIS-derived "equivalent Ay" column 
density map, to the 2MASS/NICER map in Figure UK, 
one can see that the point-to-point agreement is good, 
but far from perfect (see §3]) . Figure Q}i shows the de- 
rived color-temperature map (from ISchnee et al.l l2005h 
that is used in creating Figure [Ip. 

Note, as with the choice made for extinction maps, 
that we have not chosen to use Spitzer ma ps of thermal 
emissi on in the present study. Recently, ISchnee et al.l 
(2008) have presented a full analysis of the long wave- 
length Spitzer MIPS maps of Perseus, in concert with the 
IRAS-based data used here. The Spitzer map coverage is 
smaller than the IRAS (all-sky) coverage, and at the res- 



° Schnce ct al. (2008) have recently used the long- wavelength 
data from Spitzer to produce higher-resolution thermal-emission 
maps of column density. As we are forced to smooth to 5' resolution 
in this paper (to match the near-IR extinction maps) comparison of 
IRAS-based map s with the newer, S pitzer-based, maps is presented 
and discussed in Schnce ct al. 2008, but not here. 
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Fig. 1. — Maps in the left column show column density as traced by: Extinction (a); Dust Emission (c); and Gas Emission (c). Figure 
lb shows the small uncertainties associated with the (NICER) extinction map in Figure la. Dust color temperatures based on 60 and 
100 /im data are shown in Figure Id. Note the warm dust associated with the shell around HD278942. The opacity of 13 CO, which is 
correlated with column density, is shown in Figure If. In all panels, only data points with detections of 13 CO with S/N>= 5 are shown, 
and the resolution is 5'. The IC348 and NGC1333 regions have been excised from the data sets, because all three techniques are biased 
in cluster regions. The single pink contour surrounds the apparently heated material around HD278942, and the single light blue contour 
outlined the "overdense" area presumably created by HD278942's shell, showing the dearth of 13 CO indicated most clearly in Figure [3] 
The authors request that this figure be scaled to fill the full width of the page. 

species used in an investigation does not have a constant 
abundance relative to molecular hydrogen. In reality, any 
one a nd often all three of these are true (see lPineda et al.l 
2008) . But still, by making a prudent choice of molecular 
tracer, one can minimize these complicating conditions 

and use line maps to trace column density. 

In the past, rese archers (|Bachiller fc Cernichardll98ll 
iLanger et al.lll989| ) have argued that 13 CO, which is ex- 
cited above volume densities ~ 1000 cm~ 3 (column den- 
sities Ay ~ 1 mag in nearby molecular clouds) remains 
optically thin throughout "most" of a molecular cloud's 
volume, and also remains of relatively constant abun- 
dance, except in very cold (T < 15 K) dense (n(H2) > 
5000 cm' 3 ) regio ns where carbon-b earing species are 
heavily depleted (jCaselli et all 11999( 1. Therefore, 13 CO 



olutio n of the present study, the results of ISchnee et al.1 
(2008) show that the subtle differences between extinc- 
tions determined using 60 and 100 microns, and com- 
binations at longer wavelengths are not c ritical to our 
analy sis here (see histograms in Figure 9 of ISchnee et al.l 
l2008h . 



2.3. Molecular-Line Mapping 

The implicit assumption behind mapping column den- 
sity using molecular line emission is that by integrating 
emission over all velocities, one can trace out all of the 
molecular gas along any particular line of sight. This as- 
sumption fails when either: 1) a molecular line is only ex- 
cited under special physical conditions; 2) emission from 
a particular line becomes optically thick; and/or 3) the 
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maps are used here as the most relevant molecular line 
tracer to compare with extinction maps and IRIS-based 
dust emission maps, which are also sensitive to material 
at or above about 1 mag. 

The full 13 CO map presented here was made for 
the COMPLETE Survey at the Five College Radio 
As tronomy Ob s ervator y, and it is described in detail 
in iRidge et al.1 {2006a). The full- resolution map was 
made in on-the-fly (OTF) mode with the 32-element 
SEQUOIA array, and contains nearly 200,000 indepen- 
dent 46" pixels. As noted above, we have smoothed the 
map to 5' resolution in order to match the column den- 
sity maps described in §2.11 and §2.21 above. After this 
smoothing, the average baseline rms in each of the result- 
ing pixels is 0.027 K, and there are just over 3000 pixels 
(in all the maps in Figure [TJ. All points with signal-to- 
noise ratios (based on peak/rms temperature) below 5 in 
the smoothed map have been excised as less-than-reliable 
13 CO detections. 

One of the most common, but not necessarily most 
accurate, ways to estimate molecular gas column den- 
sity, jVfJHk), from 13 CO emission comes from: adjusting 
the velocity-integrated intensity of 13 CO, W( 13 CO) = 
J Ta( 13 CO)cIv, for telescope efficiency; assuming a uni- 
form excitation temperature; and multiplying by a con- 
version fa ctor that accounts for the ratio of H2 to 
13 CO (e.g. iBachiller &; Cernicharol 119861 and references 
therein). 

To make a more accurate calculation of iV( 13 CO) one 
can employ measured kinetic temp eratures and optica l 
depths from 12 CO observations fsee lLanger et al]|1989f ). 
Because COMPLETE includes 12 CO as well as 13 CO 
maps, we can use the observed brightness temperature 
of 12 CO to estimate the excitation temperature (assum- 
ing that 12 CO is optically thick), and then calculate a 
column density map assuming that the levels are popu- 
lated following a Boltzmann distribution. All maps and 
graphs relating to 13 CO column density distributions in 
this paper, unless expl icitly labeled as u W ( 13 COY\ are 
made in this way (see IRidge et al.l l2006afl . We also as- 
sume a constant abundance, of 3.98 x 10 5 , for H2 relative 
to 13 CO, even though some important vari ations in that 
ratio d o exist, and are studied in detail in iPineda et al.l 
(|200l . 

Figure[]J shows the opacities derived from the same set 
of calculations that yields Figure [TJ3 (the column density 
map). The opacities, many of which significantly ex- 
ceed unity, are well-correlated with the column density. 
Once t>1, any spectral line is no longer a very faithful 
column density tracer, even when we work arduously to 
correct for opacity. 

In order to convert molecular column density, N, to 
the "equivalent Ay" units used in this pap er, we use the 
procedure outlined in lPineda et al.l (|2008l ). which gives: 

A v ( 13 CO) = 4.24 x 10" 16 iV( 13 CO) + 1.67 . (1) 

To facilitat e comparison of the p resent analysis with pre- 
vious work. IPineda et all J2008) assume a constant ratio 
of reddening to extinction, Ry — Ay jE{B — V) = 3.1, 
equal to the measured ave rage value measured for the 
Galaxy (jBohlin et al.lll978l). We note, though, that Ry 
can h ave values up to ~ 6 (jDraine 2003; Goodman et al.l 
|1995|) . especially in high-density regions. Lastly, we point 
out that the coefficients in eq. U) apply to Perseus as 



a whole, even though we measure them to vary by as 
much as ±30% fro m region to region within Perseus 
(|Pineda et al.ll2008T ). We discuss the significance of these 
regional variations in §4.31 below. 

In the bottom panel of Figure [2j we show both the 
PF( 13 CO)-based and the iV( 13 CO)-based histograms of 
column density. 6 Notice that the W( 13 CO) distribu- 
tion appears to underestimate the amount of material 
at low column density (where temperature is typically 
higher) and at very high column density (where r is 
high). Even after accounting for temperature and opac- 
ity variations though, the 7V( 13 CO) distribution is still 
dangerous to interpret too literally, because at low Ay 
(<3 mag) the gas is sub-thermally excited, and at high 
extinction the optic al depth is significantly larger than 1 
(|Pineda et al.l l2008) . We quantify and discuss the uncer- 
tainties in molecular column density as compared with 
dust measurements below in §4.11 

3. DATA EDITING 

In order to restrict the comparison in this paper to re- 
liable measurements of column density in the portions 
of molecular clouds not dominated by the localized heat- 
ing and stirring caused by embedded massive stars we 
have excluded certain positions from all three data sets 
compared here. The "exclusion criteria" are as follows. 

Points not having reliable measurements in all three 
tracers are omitted. For 12 CO and 13 CO data, we use 
only positions with signal-to-noise (peak-to-rms) ratios 
greater than 15 and 5, respectively. Note, though, that 
even with infinite sensitivity there would still be a min- 
imum detectable column density for the line measure- 
ments, due to the fact that a critical density of matter is 
needed to excite CO or 13 CO collisionally. (This "mini- 
mum" density for excitation is what leads to the additive 
constant in the relationship between Ay and A( 13 CO) 
shown in eq.. {l}. We also exclude any points where the 
fitted 12 CO line-width is smaller than 80% of the 13 CO 
line-width, because this is indicative of either: 1) multi- 
ple velocity components along the line of sight captured 
by one tracer but not the other; or 2) of an unphysical 
result caused by marginal data quality. 

The IC348 and NGC1333 cluster regions are elimi- 
nated, using purely spatial cuts (see black "holes" in 
Figure QJ, from our analysis, because when material is 
heated non-uniformly from within, it is difficult to derive 
accurate column densities fr om either dust or gas emis- 
sion (see ISchnee et al.l 12006). We also remove a small 
number of other pixels with stellar densities high enough 
(in this case, > 10 stars/pixel) to imply a significant con- 
tribution from stars embedded in Perseus. This is nec- 
essary because NICER relies on "background" stars for 
extinction mapping, and will underestimate extinction if 
embedded stars are included by accident. 

As indicated on the dust temperature map in Figure 
QJ1, a heated shell surrounds the B-star HD 278942, and 
portions of it overlap with our study region, in projec- 
tion. The shell is actually located just behind the molec- 
ular clouds, and it apparently tou ches the clouds at a 
few localized points of contact (see IRidge et al.l l2006bf) . 

6 Throughout the paper, we quote column density in units of 
Ay, but there are assumptions used to reach such units, which are 
explained in J2] 
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TABLE 1 
Fit Parameters 



Tracer 


Scaling 


log A vfi 


A v ,o 


17 


A v 

Ndust 

7V( 13 CO) 
W{ 1A CO) 


336±14 
416±12 
324±27 
298±18 


0.397±0.008 
0.433±0.005 
0.40±0.02 
0.45±0.01 


1.18±0.01 
1.145±0.005 
1.19±0.02 
1.22±0.02 


0.163±0.008 
0.136±0.005 
0.17±0.02 
0.20±0.01 



In Figures Q] and El we have marked the effects of the 
shell's low-density hot interior (in pink), as well as of its 
over-dense heated rim (in blue). In the histograms of Fig- 
ure [2] and for the corresponding fits shown in Table [1] the 
points effected by the shell's interior and exterior have 
been excluded. In Figure [3l we show, for illustrative pur- 
poses only, the shell interior/exterior points in pink/blue, 
but we calculate the scatter amongst the different tracers 
excluding the shell points. 

Lastly, positions close to the border of the map are 
removed when they would be effected by inaccurate con- 
volution. 

In summary, all of the positions shown with non-null 
values in Figure [1] are shown in Figure 02 but the points 
within the "shell" regions are excluded from the his- 
tograms and fits in Figure [2j 

4. RESULTS: THE DISTRIBUTION OF COLUMN DENSITY 
IN PERSEUS 

Figure [2] compares, as histograms, the distributions of 
column density implied by the three different techniques 
(extinction, thermal emission, and molecular lines) un- 
der study here. To facilitate comparison with theories 
that predict a log- normal density distribution (see ij5.ip , 
the histograms are shown in linear-log space, where a 
log-normal appears as a Gaussian. The log-normal least- 
squares fit to the 2MASS/NICER-based extinction dis- 
tribution is repeated (in blue) as a fiducial in all plots. 
Table [1] shows the fit parameters for each column density 
distribution based on the functional form: 



Number of pixels = Scaling * e 



-(log Ay-logA V| o) /2o- 2 



It might seem surprising, especially given the non- 
gaussian shape of the 13 CO histogram, that the Ay,® 
parameter for all the column density tracers are so sim- 
ilar, however, keep in mind that this is to be expected 
because thermal emission and molecular line data have 
been calibrated with the 2MASS /NICER-based extinc- 
tion map. 

4.1. Uncertainty Based on Tracer Inter-Comparison 

Figure [3] shows direct inter-comparisons of the all three 
column density measures: extinction, thermal dust emis- 
sion, and gas emission. The overall scatter around the 
1:1 lines can likely be attributed, in large part, to vari- 
ations along the line of sight of: dust properties; dust 
temperature; and gas volume density (which effects line 
emission). We explain below that the principal source 
of variation between the two dust tracers is variation in 
temperature along the line of sight; while for comparisons 
between 13 CO and dust tracers, variations are caused by 
changes in the ratio of gas to dust along the line of sight. 

Dust emission is dependent on temperature, but ex- 
tinction is not. By applying radiative transfer cal- 




0.5 

log(Equivalent A v (mag)) 

Fig. 2. — Pixel- by-pixel column density histograms for the full 
Perseus COMPLETE data set shown in Figure Q] Note that all 
tracers arc smoothed to a common 5' resolution, and only points 
where 13 CO is reliably detected are included (see text). Each distri- 
bution shows the column density as labeled on the plot, converted 
to units of Ay, as explained in the text. In every panel, the solid 
smooth blue curve shows the Gaussian fit to the 2MASS/NICER- 
based distribution, for reference. The red smooth curve in the 
middle panel, and the green smooth curve in the lower panel, show 
Gaussian fits to the IRAS- and 3 CO-implied distributions, re- 
spectively. The grey shaded histogram in the bottom panel shows 
WT 13 CO) converted to units of Ay, but the fit shown is only for 
the curve just labeled " 13 CO," which gives total column density 
calculated using eq. JTJ. Fit parameters for all four column den- 
sity distributions shown in this Figure are given in Table [T] The 
dashed vertical line extending through all three panels shows the 
cutoff (lowest) value of column density measurable with 13 CO, 
according to eq. JTJ. The short horizontal bar centered on the 
dashed line in the top pane shows the lcr spread in distribution of 
\A V (2MASS) - Ay( 13 CO)|/A v ( 13 CO) and for the middle pane 
the same dispersion but for | Ay (IRAS) - Ay ( 13 CO)]/A v ( 13 CO)|. 
Figure l5l shows the regional breakdown of these same histograms. 

culations to nume rical models of molecular clouds, 
ISchnee et al.l (|2006f ) showed that most of the scatter in 
comparing extinction and emission based column den- 
sity measures is likely caused by line-of-sight variations 
in dust temperature. Specifically, the scatter in the mid- 
dle panel of Fig ure [3] is direct l y mod eled, and explained, 
by this effect in lSchnee et al.l (|200 6j. so we will not dis- 
cuss it further here. 

To discuss the variations between dust- and gas-based 
measures of column density, we use Figure|4l which shows 
a schematic diagram of gas and/or dust structures viewed 
from various vantage points. The general idea of a "dust- 
to-gas" ratio is fine as long as one takes care to specify 
the volume over which that ratio is intended to apply. In 
other words, if clouds like those shown in the Figure [5] 
made up a whole galaxy, then summing all the dust and 
all the gas mass and dividing would be a fine and accurate 
expression of a volume-averaged dust-to-gas ratio. In 
our case, however, where we are interested in detailed 
pixel-by-pixel comparisons of column density measured 
with either dust or gas, plane-of-the-sky variations in the 
dust-to-gas ratio matter greatly, but are difficult to treat. 

The full "dust-to-gas" ratio derivable from measure- 
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Fig. 3. — Inter-comparison of measured column density distribu- 
tions, all shown in units of Ay (see text). The horizontal dotted 
lines in the top and bottom panels show the minimum column den- 
sity measurable with 13 CO corresponding to the additive constant 
in eq. JTJ- Note that the grey scale color of the data points shows 
the "third" measure of column density not plotted on the x or y 
axis for each plot, so that any points that do not look to be part of 
a smooth grey gradient are "outliers" in that third measure. Points 
colored blue and pink correspond to the blue-outlined shell exte- 
rior and the pink-outlined shell interior in Figure [l] The pink and 
blue point are shown here for illustrative value, but they are not 
included in the histograms in Figures [T] and [5] or in any fits to col- 
umn density distributions, as explained in the text. The 45-degree 
lines are not fits: they simply show a 1:1 relationship that might 
reasonably be expected given that the IRAS-based and 13 CO-based 
measures have been calibrated to best match the 2MASS-based Ay- 
distribution overall (see text). 
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Fig. 4. — Schematic diagram showing various lines of sight 
through various conditions. A, B, and C, are views as seen from 
our vantage point on Earth. In A, we have the favorable situation 
where 13 CO is in near-LTE (shaded purple), and all of the low- 
density material around it, some of which does not emit in 13 CO 
at all (shaded grey), but still shows up in dust measures, is associ- 
ated with the cloud of interest. In B, the situation is as in A, but 
an additional region which emits in dust but not in 13 CO (either 
due to low density or due to a dearth of 13 CO) is included. In C, 
the densest gas the line of sight passes through (shown as marbled 
grey) is dense enough to excite some 13 CO, but not at a level truly 
indicative of the full column density present, because the collisions 
at this "sub-critical" density are too infrequent. In D, we see an 
"alien's view" of the same cloud, which passes through the "dust- 
only" zones but no 13 CO emitting regions, even though it crosses 
A, B, and C. 

ments of CO isotopes and extinction relies on the wave- 
length dependence of extinction and on the abundance 
of CO isotopes i n the gas. As explained in 32.31 and 
in iPineda et all (|200l . we assume in this paper that 
the wavelength dependence of extinction is constant (as 
Rv = Av/E(B — V) = 3.1), and that the abundance 
of 13 CO is fixed for all of Perseus, at the value corre- 
sponding to the coefficients in eq. ([I]). We make these 
assumptions in order to show how gas and dust maps 
would typically compare, as this approach, or an even 
less customized one, where dust-to-gas ratios and abun- 
dances from the literature are used without testing, is 
most common. In lPineda et al.l (|2008l ) we give a detailed 
analysis of how much the calibration of dust-to-gas re- 
lationships ( 13 CO abundance) varies amongst the sub- 
regions of Perseus under both the linear approximation 
(eq. ([1])) and also with a curve-of-growth approach. 

The 13 CO(1-0) transition is only excited above a cer- 
tain critical density, and the transition to "LTE" levels 
of excitation is actually a gradual one. So, in Figure [H 
we show some low-density dusty regions where 13 CO is 
not excited at all, some slightly-higher-density regions 
where it is sub-thermally excited, and some where the 
density is high enough for near-LTE conditions. Figure 
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[4] clearly shows that while we can excise points with no 
13 CO emission on the plane of the sky, we cannot excise 
dust along the line of sight that would produce no 13 CO 
emission if viewed "from the side" (e.g. by the alien in 
the Figure). 

If all lines of sight in a region were exactly the same 
as line-of-sight "A" in Figure 01 then: 1) an equation for 
converting gas to dust column density like eq. could 
be valid; and 2) there would be no scatter in comparing 
dust and gas measures of column density. However, in a 
real cloud like Perseus, some lines of sight pass through 
additional material that contains dust but does not emit 
in 13 CO, as shown for line of sight "B" in Figured] That 
additional material can be background or foreground to 
the cloud. Or, it can be material within the cloud that 
is depleted or deficient in 13 CO, for a variety of reasons. 

The 13 CO column densities given in this paper 
are calibrated using a single "typical" conversion fac- 
tor in eq. (JTJ) , for all of Perseus. This calibra- 
tion makes the implicit assumption that the mix- 
ture of foreground/background/depleted/deficient mate- 
rial along lines of sight within Perseus does not vary too 
much, which is not strictly true-and the scatter seen in 
Figures [2] and [3] is caused primarily by this assumption. 
Note that since we know that the shell region (highlighted 
using blue and pink in Figures [T] and [3 is unusual in both 
its temperature and for its apparent 1 CO deficiency (see 
i j4.21 below), it has been excised from all fits (including 
the overall calibration) and histograms (Figures [T] and [5]) 
in this paper. 

Equation (TTJ) has an additive factor corresponding to 
the typical minimum column density at which 13 CO is 
detected. This minimum is shown with dotted lines in 
Figures [2] and E] at logA v = logl.67 = 0.22. Just as the 
slope of the gas-to-dust calibration does not universally 
apply to all lines of sight, neither will this minimum. 
As explained above, we have excised all positions with 
13 CO levels below the minimum set by eq. (fTJ), but the 
remaining positions will sometimes give dust-derived ex- 
tinctions below the gas-derived minimum. The biggest 
effect causing "leakage," beyond the low-density cutoff, 
which is partially responsible for the low-density tails 
seen in Figures [2] and [5j is simply the scatter about the 
1:1 relationships shown in Figure [3l This scatter is char- 
acterized by horizontal error bars in Figures [5] and [5l the 
calculation of which is explained below. 

A more subtle effect causing "leakage" below the 13 CO 
low column cutoff is also present, and it is illustrated 
in line of sight "C" in Figure [4] which shows a line of 
sight that passes through a "sub-thermally excited" re- 
gion. In cases like "C," 13 CO will be detectable, but 
the density in the region producing the emission is below 
or barely at the critical density, such that one cannot 
reliably convert 13 CO line intensi ty to a column den - 
sity using eq. (TTJ). As shown in iPineda et al.l (2008), 
when one uses a more appropriate curve-of-growth fit, 
rather than a linear approximation, 13 CO integrated in- 
tensity typically rises more steeply than a linear fit to 
the same data would imply. As a result, a linear con- 
version (eq. <(T])) applied to the observed 13 CO intensity 
in low-density (sub-thermally-excited) regions will often 
overestimate the column density. Thus, some of the posi- 
tions creating the "low column" tails shown in Figures [2] 
and[5]are likely to actually be at the low column densities 



that dust measures: their column density is just overes- 
timated by (sub-thermally-excited) 13 CO. Note, further, 
that the seemingly odd pile-up of 13 CO at densities just 
above the cutoff in Figure [2] is similarly caused by the 
linear approximation's inability to properly treat sub- 
thermally-excited (relatively low density) gas. 

To empirically estimate the scatter associated with 
all manner of variations in dust and gas properties, 
we calculate the standard deviation of the distribu- 
tions of the normalized differences between one col- 
umn density measure and another at each point in 
the maps. We find, for the full Perseus map, that 
that |A y (IRAS)-Ay(2MASS)|/A y (2MASS) has 1-sigma 
width (standard deviation) of 26% and \A v ( 13 CO)- 
^y(2MASS)|/^ y (2MASS) has la width 24%. It is clear 
from Figure [3] that these standard deviations will be 
larger than the la uncertainty in a linear fit, because 
the 1:1 line (implicitly assumed in the standard devia- 
tion calculation) is not a perfect representation of the 
real relationships between all the extinction measures. 

In Figures [5] and we show horizontal bars of vary- 
ing lengths centered at the threshold value of log Ay — 
0.22 in order to demonstrate that the fluctuations one 
sees about the 1:1 lines in Figure [3] are large enough 
to cause the low-column-density tails seen in Figures 
[2] and [5l The specific length for each horizontal 
bar shown is the standard deviation of |^4y(IRAS)- 
A v ( 13 CO)\/A v ( 13 CO) for IRAS-based panels of Figures 
HandEl and \A v {2MASS)-A v { 13 CO)\/A v { 13 GO) for the 
2MASS/NICER-based panels. We discuss the implica- 
tions of dust and gas probing slightly different regions 
along the line of sight further in £15-31 below. 

4.2. Paucity of 13 CO in the Shell Around ED 278942 

The overabundance of blue points in the lower-right 
portion of the two parts of Figure [3] that involve 13 CO 
data is likely due to a dearth of 13 CO in the shell around 
HD 278942. On the plane of the sky, nearly all of the 
points that lie below the overall 1:1 trend are associ- 
ated with the shell (see blue contour in Figure [1]). The 
effect appears most strongly when comparing 13 CO to 
dust emission (Fig. [3l top), but we base our statements 
here on the comparison of 13 CO and extinction (Fig. [3l 
bottom), because the shell's effects also biases the dust 
emission column density measurements in that region (as 
shown in the middle panel of Figure [3]). We suspect that 
either the shell structure is young, and molecules, such 
as 13 CO, have not had time to form, or th at the ener- 
getic r adiation associated with the shell (see lRidge et al.1 
2006b) has destroyed pre-existing 13 CO. In other words, 
the column density along lines of sight through the shell 
(outlined in blue in Figure[T]) has a molecular component, 
associated with the (rest of the) Perseus star-forming re- 
gion, plus a component associated with the shell which 
contains a lower fraction of molecular gas. 

4.3. Regional Variations 

Perseus itself is a large complex of molecular clouds. 
A division of it into a set of morphologically-distinct 
sub-regions can be guided by a plo t of dust temperatur e 
versus gas velocity, as explained in IPineda et al.l ([2008) . 
Histograms of column density analogous to those shown 
in Figure [5] are shown in Figure [5] for the same sub- 
regions (B5, IC348, Shell, Bl, NGC1333, and Westend) 
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Fig. 5. — Regional variations in column density distributions. Each panel shows the same assortment of measured distributions as in 
Figure [2] only without VF( 13 CO). As explained in more detail in Figure [2]s caption and in the text, the vertical dotted line shows the 
minimum Ay measurable with 13 CO and the short horizontal bars show indicate the characteristic scatter about the relationship between 
each dust tracer and 13 CO. Parameters for Gaussian fits to normalized versions of these distributions are given in Table[2] The short vertical 
lines hanging from the top axis indicate the mean value of extinction for each distribution (by color), and the means of the normalized 
extinction are given in Table[2] Note that we use "frequency" on the y-axis in these plots, rather than "number" as in Figure[2] to facilitate 
inter-comparison of the various regions, not all of which include the same total number of positions (see Table [2} • 



discussed in iPineda et al.l ([2008). Significant variation 
from region to region is clearly present. Typically, each 
sub-region has of order 500 points (out of ~ 3500 for all 
of Perseus, see Table [5]) and has a maximum linear extent 
of ~ 5 pc rather than the ~ 25 pc full length of Perseus, 
making each sub-region a significantly smaller statistical 
and physical sample than is all of Perseus. 

The shapes of the distributions vary significantly from 
sub-region to sub-region, and most also show significant 
disagreement amongst the three column density tracers 
used here. One key point is that some of the tracer-to- 
tracer disagreement is caused by (purposely) using only 
a single CO abundance, and a single form of eq. |1| in 
creating Figure O Had we cu stomized the calibr ation 
region- by-region, as is done in IPineda et alj (|2008l ). we 
could improve the tracer-to-tracer agreement some, but 
then we would not be representing legitimate sub-sets of 
the same data shown in Figure Q] through [3l Overall, 
13 CO as a tracer seems most capricious, and we suspect 
that this is due to significant variations in physical condi- 
tions (e.g. temperature, radiation field, region age) other 
than (column) density that effect its abundance and ex- 
citation. 

A subtler point concerning variations in distribution 
shape, both regionally, and tracer-to-tracer, concerns 
real variations in dust-to-gas ratio, in gas properties 
(e.g. 13 CO abundance), and/or in dust properties (e.g. 
i?y). We know that these properties vary (cf. H2.3I and 



IPineda et all 12008) even on the many-pc scales charac- 
teristic of the "regions" we discuss here. So, they likely 
also vary on smaller scales, meaning that subtle changes 
in the shape of the distributions seen in Figure [5] and 
O would be apparent if we could account for these vari- 
ations. In th e ir pap er on "Can We Trust the Dust?," 
iPadoan et al.l (|2006t see also t j5.3l below) suggest that 
very small-scale variations (< 0.1 pc) variations in the 
dust-to-gas ratio may in fact exist. So, in the future, 
when we can carry out inter-comparisons of column den- 
sity tracers on even smaller scales than those we con- 
sider in this paper, it will be interesting-and potentially 
important-to quantify how much the dust-to-gas ratio, 
as well as intrinsic gas and dust properties, change at 
various scales, and under various conditions. 

If we study only the NICER-bascd histograms in Fig- 
ure [5l which are not effected by calibration choices here, 
there are still very significant variations from region to 
region, and one is left wondering how sample (region) size 
effects the shape of one of these histograms. We consider 
this question, and others, in the context of numerical 
simulations, below. 

5. IMPLICATIONS 

5.1. Comparisons with Numerical Simulations 

In this section, we consider how numerical simulations 
offer insight into the shape of column density distribu- 
tions. 
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TABLE 2 
Regional Variations 



Name Mean Ay Mean In Ay Mean In x Sigma In a; 



All, Number of Points= 2892 



NICER 3.080 1.024 -0.101 0.433 

IRAS 3.032 1.052 -0.057 0.329 

N( 13 CO) 3.243 1.076 -0.101 0.427 

W( 13 CO) 3.215 1.083 -0.084 0.404 



B5, Number of Points= 344 



NICER 2.734 0.940 -0.065 0.349 

IRAS 2.692 0.968 -0.022 0.213 

N( 13 CO) 2.850 0.975 -0.073 0.368 

W( 13 CO) 2.733 0.941 -0.065 0.353 



IC348, Number of Points= 552 



NICER 3.872 1.230 -0.124 0.500 

IRAS 3.826 1.279 -0.063 0.359 

N( 13 CO) 3.487 1.105 -0.144 0.525 

W( 13 CO) 3.293 1.068 -0.123 0.496 



Shell, Number of Points= 176 



NICER 3.318 1.076 -0.123 0.503 

IRAS 3.523 1.194 -0.065 0.373 

N( 13 CO) 3.250 1.053 -0.125 0.484 

W( 13 CO) 3.193 1.053 -0.108 0.458 



Bl, Number of Points= 631 



NICER 3.268 1.075 -0.109 0.448 

IRAS 3.057 1.068 -0.050 0.304 

N( 13 CO) 3.351 1.086 -0.123 0.477 

W( 13 CO) 3.233 1.069 -0.104 0.450 



NGC1333, Number of Points= 642 



NICER 2.681 0.915 -0.071 0.358 

IRAS 2.813 0.994 -0.041 0.271 

N( 13 CO) 3.445 1.147 -0.090 0.396 

W( 13 CO) 3.518 1.190 -0.068 0.356 



Westend, Number of Points= 547 



NICER 2.674 0.919 -0.065 0.353 

IRAS 2.511 0.879 -0.041 0.285 

N( 13 CO) 2.879 1.021 -0.036 0.259 

W( 13 CO) 3.071 1.089 -0.032 0.250 



IVazquez~ Scmadcnj| (|1994| ) shows that for highly su- 
personic flows where gravitational and magnetic forces 
become negligible, the gas has a pressureless behavior. 
Under these conditions, the hydrodynamic equations be- 
come scale invariant, i.e., motions at all length and den- 
sity scales obey the same equations. As a result, the 
probability density functio n of the volume den sity, n, is 
expected to be log-normal. lOstriker et al.l (|2001l ) demon- 
strate that for essentially isotropic flows, the same kind 
of log-normal distribution results for either volume or for 
column density. These general predictions of turbulence 
theory inspire the log-normal fits shown in Figure O We 
do not suggest that any of the column density distribu- 
tions we find here are necessarily best-fit by log-normals- 
we simply offer log-normal fits as a relevant comparison. 



Figures [T] and [5] straightforwardly show number (Fig- 
ure [I) and frequency (Figure 0) distributions of logA^. 
If, however, one wants to emphasize density fluctuations 
about a mean, it makes sense to normalize the distribu- 
tions by a mean column density. Most numerical sim- 
ulations, because they are often scale invariant, analyze 
this kind of "normalized column density," so to facilitate 
comparison, we define 

x = N/N (3) 

where N is the mean value of column density in any 
map. The values of such means are given in Tableland 
shown as vertical long ticks in Figure [5] for all the data 
sets considered here. 

We define the standard deviation of the distribution of 
In a; as o~i nx , which would equal the la standard deviation 
of a Gaussian fit to the distribution of In x in the case of a 
pure log-normal column density distribution. The values 
listed for o~\ n x in Table [5] are determined purely from the 
statistics of the distribution of \nx values, though, and 
not from any kind of (e.g. Gaussian) fit. 

Some simulations have suggested that there may be 
measurable relationships amongst the mean and disper- 
sion of the distribution of In x, Mach number, mean mag- 
netic field strength, and the ratio of forcin g scale to cloud 
scale. For example, iPadoan et al.l (|1997bT ) suggest that 

«r ln;c 2 = m(l + 7W 2 /? 2 ) (4) 

where (3 is a constant of order 0.5 and M. is the sonic 
Mach number of the gas. However, if a m agnetic field is 
presen t, both IPadoan et al.l (|1997bf ) and lOstriker et al.l 
(120011) find that this relationship changes. In particular, 
lOstriker et al.1 (|2001h find a secular trend in density con- 
trast that depends on the fast magetosonic Mach num- 
ber, Mf, which in turn depends on the sound speed and 
the Alfven speed. Presently, region-by-region detailed 
measurements of field strength, needed to measure the 
Alfven and sound speed independently, are not available 
for any set of regions as large as those we study here, so it 
is hard to test this relationship directly with observations 
right now. However, a concerted effort to test predictions 
that rely oiiMf, perhaps using the Chadrasekhar- Fermi 
method to estimate field strengths over large areas, could 
and should be undertaken in the near future. 

Even if we could know field strengths though, it is still 
not the case, according to simulations, that a single real- 
ization of a turbulent flow can be inverted to give basic 
physical parameters, such as the exact power spectrum 
of d ensity fluctuations in the flow. Several resea rchers 
(e.g. IVazquez-Semadenilll99l lOstriker et ail 1200 ll ) have 
shown that sub-samples of a single flow and/or multi- 
ple realizations of the same physical conditions will give 
noticeably different (normalized) column density distri- 
butions, similar to what is seen in Figure There are 
a large number of modes present simultaneously in any 
simulation so one sample alone is unlikely to sample them 
all in their "average" mixture. 

So, given that simulations demonstrate that a large(r) 
sample of realizations would be needed to causally re- 
late basic physical parameters to column density distri- 
butions, what can we learn from just the data we present 
in this paper about the physics of turbulent flows in star 
forming regions? 

We can consider our set of "regional" sub-samples as 
multiple realizations of the same "experiment" we can 
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Fig. 6. — Comparison of the fit parameters listed in Table[2]with 
each other, and with the 1:1 line predicted, for a log-normal, by 
eq. 

call the "Perseus Molecular Cloud." In one sub-sample, 
the "Shell" region, where the gas is dominated by an ob- 
vious driver not present in the other regions, we might 
dismiss any outlier-like behavior as due to "unusual" 
forcing, and there is in fact a skewing of its column 
density distribution toward lower values, not seen in the 
other regions (see Figure [5]). 

If the regional density distributions are truly drawn 
from an underlying distribution which is inherently log- 
normal, as predicted by many simulations, then 

2 

m^=-^i ( 5 ) 
2 v ; 

simply because the mean (first moment) of a log-normal 
distribution of a quantity (e.g. x) that is normalized by 
its own mean should be unity. Thus, Figure [6l which 
shows — c^ lx /2 as a function of the mean lnx for each 
region, appears to indicate that nearly all of the distribu- 
tions we study in this paper are close to consistency with 
being drawn from a log-normal. Note, however, that dis- 
tributions similar in shape to a log-normal will also give 
a relationship very similar to eq. ©, so the fact that 
the points for our regions lie so close to the line should 
only be taken to mean that the distributions are close to 
consistent log-normal, and not exactly log-normal. The 
main purpose of Figure [5] is to show all the values shown 
in Table [2] together, in a way that facilitates testing hy- 
potheses relying on making comparisons. 

If a relationship similar to eq. were satisfied, then 
moving up and right along the "perfect log-normal" line 
in Figure O would mean higher Mach number in the hy- 
drodynamic case, or possib ly higher fast magn etosonic 
Mach number (as shown bv lOstriker et al.1l200lL in their 
Figure 4). Note that M will rise if the sound speed drops, 
and M. f will rise either if the magnetic field drops or if 
the sound speed drops, for a given density of material. 
We have tested whether a sonic Mach number determined 



fr om observed line wid ths and (dust) temperatures listed 
m iPineda et all ([20081 ) increases along the 1:1 line in Fig- 
ure and it does not. 

Making comparisons based on Mach numbers calcu- 
lated using only dust (not gas) temperature and line- 
of-sight velocity dispersion limits direct observation- 
simulation comparison in two important ways. First, we 
simp ly cannot calculate Mf (as is used in lOstriker et al.1 
|2001| ). which is likely more relevant than M in com- 
parisons, because we do not have field strength mea- 
surements. Second, in observations, the observed col- 
umn distribution arises from plane-of-sky compressions, 
whereas the line widths used to compute Mach number 
come from the line-of-sight motions only. We consider 
how one might get around these limitations in the future 
in H5.4I below. 

One important point to take away from observation- 
simulation comparisons concerns the sometimes- 
overlooked distinction between volume density and 
column density distributions. Several simulators have 
pointed out that the column density distributions 
we show here will only look the same as the volume 
density distributions, which are usually predicted to be 
log-normal, w hen turbulence is do minated by large-scale 
motions (e.g. lOstriker et al.l [20011 ) . The fact that the 
distributions look as close to log-normal as they do 
suggests that turbulence in Perseus is primarily driven 
by motions on scales nearly the size of the whole cloud, 
or larger. 

5.2. How common are log-normal-like density 
distributions ? 

A similar relationship to the one shown in Figure O 
between the mean and dispersion in column density, 
has been used to explain and interpret the results of 
lLada et ail (|1994f ) , who we believe were the first to find a 
relation between the standard deviation i n measured ex- 
tinctio n and the extinction value itself. [Thorav al et al.l 
(1997) pointed out that a relation like eq. ^ can be 
caused by small scale structure below the re solution limit 
of an extinction map. iPadoan ~et al. (1997a) performed a 
reanalysis of Lada et al.'s IC5146 observations, and find 
that the a a v —Ay relation is consistent with a log-normal 
column density distri bution producing unr esolved den- 
sity structure. Later, lOstriker et al.l (|2001f) us ed higher 
resolution observation in IC5146 bv lLada et alj (|1999f ) to 
compare the cumulative column density distribution with 
their own simulations finding that they are very similar 
in shape. 

Other extinction studies might provide notable excep- 
tions to the "log-normal" trend. For example, the Pipe 
Nebula's column-density distribution (determined using 
the 2MASS /NICER method) displays a c omplex shape 
with multiple peaks (|Lombardi et al.ll2006| ). Lombardi et 
al. attribute these peaks to the possibility of background 
clouds observed in projection with (and thus in addition 
to) the Pipe at low Gal actic latitudes. Rece nt C 18 ob- 
servations in the Pipe (|Muench et al.l [20071 ) have shown 

7 Note that we have excised regions from "Perseus" that are 
likely to contain, and be effected by, the most local sources of tur- 
bulence: the densest portions of the NGC1333 and IC348 young 
clusters have been excluded from our analysis, so the stellar wind 
and outflow drivers they contain are not included in the area ana- 
lyzed here. 
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that there are two distinct velocity components there, 
and so it is possible that if those distributions could be 
separated, each would have a more log-normal-like col- 
umn density distribution. 

In general, it is important to consider what boundaries, 
in either spatial or velocity dimensions, are placed upon 
regions when assessing how common log-normal distri- 
butions might be. Overlapping "clouds" along the line 
of sight may appear as one structure on the plane of the 
sky, and, without kinematic or distance information, it 
is hard to separate these structures. And, even if one 
can separate populations along a line of sight, the mean- 
ing of the region one considers on the plane of the sky 
needs to be understood as physically "small" or "large" 
in the context of a particular question, such as how "log- 
normal" a distribution might be expected to look. 

5.3. Which Tracer Tells us What, how Reliably? 

Rece ntly, in a paper en titled "Can We Trust the 
Dust?". |Padoan et al.l (|2006l ) have raised questions about 
the reliability of NIR extinction mapping, based on their 
finding that the power spectrum of a NIR extinction map 
of the Taurus region is significantly shallower than that 
derived from a 13 CO map of the same region. Through 
detailed modelling, they rule out depletion of CO, and/or 
CO formation timescales as the cause(s) of this discrep- 
ancy. By producing synthetic NIR-extinction observa- 
tions of simulations, they show that the power spectrum 
derived from observations should trace the power spec- 
trum of the actual spatial distribution of dust if the ex- 
ti nction is prop o rtiona l to dust column density. 

iPadoan et all ([2006) argue that a discrepancy between 
the power spectra of the dust and the gas might be ex- 
pected in the case of transonic, nonmagnetized turbu- 
lence. Using 13 CO and NIR data as we do here, but 
for Taurus, they interpret a shallower power spectrum 
in dust fluctuations than in gas fluctuations as evidence 
for intrinsic spatial fluctuations in the dust-to-gas ratio, 
wit h amplitude in c reasin g toward smaller scales. Strictly, 
the IPadoan et al. (2006) work is only applicable at very 
small scales (their estimate is < 0.1 pc). Thus, the ap- 
parent disagreement between our results, which show the 
gas to be less reliable, and theirs which question the ve- 
racity of dust as a tracer of column density, may be only 
a question of scale. 

On the other hand, as shown bv lPineda et al.l (|2008f ) in 
Perseus, the 13 CO-dervied column density estimates are 
adversely effected by optical depth even at Ay ~ 4 mag, 
and variation in the 13 CO abundance with respect to H2 
is found between regions in the same molecular cloud. 
These caveats in the interpretation of the 1 3 CO emission 
maps are not fully taken into account in the IPadoan et all 
(2006) analysis (which assumes 13 CO to trace density 
faithfully out to Ay ~ 10 mag), and may also explain 
the apparent disagreements-both between gas and dust, 
and between our results and "Can We Trust the Dust?." 

We have clearly demonstrated, in every figure in this 
paper, that gas-based and dust-based measures of col- 
umn density almost never agree perfectly in detail. So, 
instead of "Can We Trust the Dust?" , one might ask in- 
stead "Should we Sass the Gas?". As we discussed in 
^A.ll dust and gas measures are virtually never tracing 
exactly the same portions of any particular line of sight, 
so some amount of "disagreement" when inter-comparing 



tracers has to be expected. It also has to be recognized, 
though, that any particular molecular transition can only 
trace a (sometimes very) limited range of volume densi- 
ties, which is effectively bounded at bottom by the crit- 
ical density and at top by depletion and/or opacity. Be- 
cause it is hard to know the real density range to which a 
line is sensitive, and to model how linear the conversion 
from line flux to density is likely to be, gas really does 
deserve some degree of sassing if one wants to tease its 
ability to trace column density cleanly. 

Dust, however, is not perfect either. For extinction- 
based measures, variations in the compositional and size 
distribution of grains can effect the reddening law. And, 
for emission-based measures, variations in grain temper- 
ature along the line of sight are impossi ble to correct 
for an d so impose unavoidable uncertainty (|Schnee et al.l 
l2006h . 

The principle advantage of Dust over Gas is dynamic 
range. No single observation of a gas tracer can sample 
more than about a factor of a about ten in column den- 
sity, but, dust-based measurements can span a dynamic 
range of 50 or (much) more. For example, in the famous 
B68 Globule, obs ervations of extin ction span the range 
- to 27 mag (jAlves et al.l 120011 ) while C 18 ranges 
from 0.2 to 0.8 K km s" 1 and then depletes, and N 2 H+, 
which only is detected abo ve about Ay of 10 , spans the 
range 0.3 to 2.1 K km s _1 (jBergin et al.ll2002l Figure 1). 
In TMC-1C, thermal dust emission allows for a dynamic 
range of nearly 25 in Ay (from 4 to 90 mag, iSchnee et al.l 
l2007bf) . but none of the many lin es mapped there ha s a 
dynamic range of more than ten (|Schnee et al.ll20"07af ). 

Extinction measures will always saturate at some high 
value of extinction where background sources cannot be 
detected, and emission is often limited in utility at very 
low extinctions, due to low flux levels and/or observa- 
tional strategies (e.g. chopping) that make mapping ex- 
tended emission difficult. In cases where one can can 
cross-calibrate extinction and emission using maps of re- 
gions where they should give the same information, one 
can stretch the dynamic range of "dust" to its fullest by 
using extinction as the sole probe at low column densi- 
ties, emission alone at the higher ones, and both extinc- 
tion and emission in the intermediate regime. 

In the maps of Perseus in this paper, we have excised 
points with no 13 CO emission, but our dust measures 
still can trace low-density material along those lines of 
sight that do emit in 13 CO (se e ij4.ll and Figure H]) . In 
Figure 6 of ISchnee et all ((2005) we show the dust-based 
column density distributions for all of Perseus, without 
limiting ourselves to 13 CO-emitting regions. Those dis- 
tributions look similar to the ones in Figure except 
that they show a longer high-density tail, caused in large 
part by not excluding all of the high-density cluster re- 
gions NGC1333 and IC348. It is very important to realize 
that our excision of plane-of-the sky points with no 13 CO 
here does not exclude low-column regions traced by dust 
along the line of sight, which is why the distributions in 
Figures fj] and [5] look more symmetric (not cut-off at the 
low end) for the dust measures than for 13 CO. 

The only advantage of Gas over Dust when mapping 
column density is the ability it offers to kinematically 
separate regions along the line of sight. Otherwise, Dust 
wins. 
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5.4. How to Map Column Density Distributions in the 

Future ? 

From the discussion above, we may sound ready to rec- 
ommend dust as the single best tracer of column density 
in interstellar space-but we are not. We cannot forget 
that the total (gas + dust) column density in any region 
can only be as accurate as the gas-to-dust conversion 
factor used. And, since all interstellar regions contain 
much more gas than dust, a small error in a dust-dervied 
column will be compounded when converted to total col- 
umn. 

Therefore, we recommend a "holistic" approach to 
measuring column density. We need to take account 
of the level of uncertainty inherent in each conversion 
factor (e.g. abundance ratio, reddening-to-extinction ra- 
tio, dust opacity, dust-to-gas ratio) and assumption (e.g. 
about line of sight structure) we use in a column den- 
sity calculation, and then we need to choose the best 
solution for each particular question. For example, in a 
circumstance where we are relatively confident that the 
dust-to-gas ratio is not varying, and that line-of-sight 
blending is not causing confusion, then extinction maps 
will provide the best handle on the kind of relative col- 
umn density measurements needed to measure key dis- 
tribution functions, such as "clump" mass functions. On 
the other hand, in a region where the dust-to-gas ratio is 
very uncertain (e.g. near an HII region) using molecular 
line ratios to measure volume density and opacity and 
converting to column density may be superior. In other 
words, we recommend that more care than has been typ- 
ical in the past be given to making the most appropriate 
choice in any particular study. 

Lastly, several recent studies have shown that "observ- 
ing" simulations using radiative transfer and/or chem- 
istry codes, and synthetic telescopes that mimic bi- 
ases imposed by real ones (an approach we have called 
"Taste- Testing" ) can sometimes uncover hidden limita- 
tions and biases a s sociated with certa in techniques (e.g. 
ISchnee et aHl200d IPadoan et alJl2006f >. Thus, if a realis- 
tic simulation can be observed synthetically with various 
column density probes, it can provide a good guide to 
which kind(s) of column density measures might be least 
biased under relevant conditions. 

5.5. Conclusions 

Careful re-calibration and inter-comparison of extinc- 
tion, thermal emission, and molecular emission maps of 
Perseus has allowed us to conclude that: 

• The column density distribution of material in the 
full Perseus star- forming region, with 1 < Ay < 12 
mag, is roughly log-normal, when it is not directly 
effected by embedded clusters or young stellar out- 
lows (bipolar or spherical). 

• Dust is superior to molecular lines for tracing out 
the "full" mass distribution over the range of ex- 
tinction studied, because it does not require a 
threshold density to "excite" and it does not die 
out at high densities due to high opacity or chem- 
ical depletion, the way 13 CO does. 

• The dearth of molecular gas ( 13 CO) in the region 
corresponding the shell created by the B-star HD 
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278942 suggests that either CO has not yet formed 
in this young structure, and/or that existing molec- 
ular gas has been dissociated by the shell's interac- 
tion with the cloud. 

• When Perseus is dissected into smaller "sub- 
regions," the column density distributions become 
more ragged, as is predicted by simulations for 
samples that are statistically small. However, we 
find that the sub-regions distributions are still log- 
normal-like, in that the relationship between their 
normalized means and their dispersions follows a 
trend consistent with log-normal distribution. 

• In comparing observations of column density (or 
mass) distributions with each other, and/or with 
simulations, it is perhaps more important than has 
been previously appreciated to account for the ef- 
fects of biases due to dust temperature variations, 
abundance variations, opacity effects, and observ- 
ing strategies. 

We recommend, for the near-term future, the assem- 
blage of an ensemble of maps of column density, along 
with measurements of Mach number and magnetic field 
strength, in order to assess the turbulent properties of 
molecular clouds and star-forming regions more gener- 
ally. This large sample is needed to allow for legitimate 
comparisons between regions, and of observations with 
simulations, because any one observation of even the 
same turbulent flow is not enough to characterize the 
flow's statistical nature. With such a large sample, we 
could carry out much more discriminating comparative 
analyses than the kind represented by Figure El which fo- 
cuses only on testing whether the conditions relating the 
mean (first moment) of a distribution to its width (sec- 
ond moment) are correct. We could, for example, begin 
to investigate the skewess of these distributions, and to 
investigate alternative functional forms, which may in 
fact not be exactly log-normal. With this larger sam- 
ple, we could also study the effects of star-formation as 
a driver of the turbulence, and again compare with sim- 
ulations of this process. 

The combination of: 1) extensive recent improvements 
in extinction mapping made possible by large-scale near- 
infrared surveys; 2) the advent of huge molecular line 
surveys of relatively high-density tracers; and 3) in- 
creases in polarization mapping speed (which leads to 
Chandrasekhar- Fermi-based field estimates) , should very 
soon allow for studies large enough to test predictive the- 
ories and simulations of molecular cloud topology using 
unprecedentedly large observational statistical samples, 
and we look forward to it. 

This paper was originally inspired by a workshop at 
the Aspen Center for Physics in the Summer of 2004 on 
"Star Formation in Galaxies," where it seemed that the 
assembled audience of experts could not agree on the 
least biased way to measure the "initial conditions" for 
stars to form from molecular gas. A conversation with 
Eve Ostriker at that meeting, about how observations 
and simulations of star-forming molecular gas might best 
be compared, was particularly important. The quest to 
offer the most bias- and error-free column density distri- 
butions we could publish based on the COMPLETE data 
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in this paper took nearly four y ears, and it spawned sev- 
eral other papers by our group (ISchnee et al.ll200(j 120081 : 
Rid ge et all2006at Pineda et al Jl2008HFoster et al.ll2008h 
The "2007" version of these distributions and their im- 
plications were discussed intensively at the KITP Santa 
Barbara Workshop on "Star Formation Near and Far," 
and we deeply thank Eve Ostriker, Paolo Padoan and 
Enrique Vazquez-Semadeni for their comments both at 
and since that meeting. We thank Joao Alves, Michelle 
Borkin, Paola Casern, Jonathan Foster, Jens Kauffmann, 
Di Li, Marco Lombardi, and Naomi Ridge for their im- 
portant contributions to the data and results presented 
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